Electric field modulation of topological order in thin film semiconductors 
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We propose a method that can consecutively modulate the topological orders or the number of 
helical edge states in ultrathin film semiconductors without a magnetic field. By applying a staggered 
periodic potential, the system undergoes a transition from a topological trivial insulating state into 
a non-trivial one with helical edge states emerging in the band gap. Further study demonstrates 
that the number of helical edge state can be modulated by the amplitude and the geometry of the 
electric potential in a step-wise fashion, which is analogous to tuning the integer quantum Hall 
conductance by a megntic field. We address the feasibility of experimental measurement of this 
topological transition. 
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Edge and surface physics arising from the topological 
insulators have been a major research focus in the con- 
densed matter physics recently. 0, i, 0, fl 0, 0, 0, 0, H 
Usually spin-orbit coupling in these materials is so strong 
that the band gap between the conduction and valence 
bands can be inverted. In this case edge states or sur- 
face states emerge in the band gap, protected by time- 
reversal symmetry. The Z 2 invariant is established to 
govern the topological properties of both two- and three- 
dimensional (3D) topological insulators. If], 11 1 Search- 
ing for this topological state of quantum matter has been 
focused on specific materials such as HgTe/CdTe quan- 
tum wells Bismuth alloys jf| 0, 0, 0] and transition 
metal oxide Na2ir03.[12|] It was observed that there are 
five surface states in Bii-^SbxQ and single Dirac cone in 
Bi 2 Se3 0] and Bi 2 Te 3 [9j] while there is only a pair of helical 
edge states in HgTe/CdTe quantum wells. 0, [l3j] Up to 
now there is neither experimental report nor theoretical 
proposal to tune the topological number in these system 
in a controllable way. It would be both theoretically and 
experimentally interesting if the topological number or 
the number of the helical edge states or surface states 
can be modulated like the quantum Hall conductance in 
a magnetic field. 

In this paper, we propose a feasible approach to mod- 
ulate consecutively the pair number of helical edge states 
by an electric means. To be specific, we consider a quasi 
two dimensional semiconductor thin film or quantum well 
with strong spin-orbit coupling, which can be either topo- 
logically trivial or non-trivial. We apply a lateral surface 
superlattice (SSL) that creates a periodic potential on 
the film. It is found that the system changes from a triv- 
ial insulator into a quantum spin Hall (QSH) insulator 
with helical edge states emerging at the sample edges. 
This phenomenon can be understood intuitively as the 
periodic potential splits the band structure into multi- 
mini-bands, inverts the band gap and changes the Z 2 
order of the bulk. However, unlike the QSH insulator, 
the band inversion story doesn't stop here. By further 



increasing the gate geometry or potential magnitude, the 
pair number of helical edge states goes up step wise while 
the bulk band gap remains finite. In this manner we cre- 
ate a topological insulator whose edge states are tunable 
by purely electrical means. 

We start with an effective four-band Hamiltonian with 
the time-reversal invariance, 



Hn = 



h + (—id x , -id y ) 






h-(-id x , -idy) 



(1) 



where h± = +D(d% + d?)+A(-id x a y +id y a x )±[A/2 



B(d% + dy)]a z and a a are the Pauli matrices. The model 
was first introduced for the HgTe/CdTe quantum well 
for QSH effect by Bernevig, Hughes and Zhang, Q and 
recently was derived for an ultra-thin film of 3D topo- 
logical insulator. [15( The two cases have different basis 
although the forms are identical. h± consists of the in- 
variants in the irreducible representation D 1 / 2 of SU(2) 
group. Hi The whole model Hq keeps the time reversal 
invariance. h + and h_ are the "time" reversal counter- 
parts under the operation O = —io y K where K is the 
complex conjugate operator, h + = Q h-Q. The block 
diagonalized form in Hq allows us to study h+ and h— 
separately and then put together to gain the physics for 
Ho. Additional terms such as those of bulk or structure 
inversion asymmetry can couple them together. Here we 
first assume that these effects are weak and negligible. 

The idea of a lateral SSL dates back to the 1970s. Var- 
ious alternative techniques are known capable of creat- 
ing SSLs.[3] It is known that periodic potential created 
by SSLs induces Bloch mini-bands and mini-gaps on a 
2DEG, which results in interesting transport behaviors. 
Consider this model on a SSL. To be specific and with- 
out loss of generality, we introduce a square- wave-shaped 
periodic potential, 

v(v) = ( Vo > °^y< d / 2 (2 ) 

1 -V , -d/2 ^y <0 [Z) 
and V(y + d) — V(y). We write the Fourier series for 
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FIG. 1: Evolution of energy bands in a modulated periodic 
potential, (a) The band in the absence of periodic potential; 
(b) Mini-bands formed in a periodic potential; (c) Shifting 
of the minibands leads to the band invertion and band anti- 
crossing. 

V(y) as V{y) = V n e m ® v with the reciprocal vec- 

n 

tor Q = 2n/d and the Fourier transform component 
V n = 2iVo/nir Jn = ±1, ±3, ■■■.)■ According to the 
Bloch theorem, [17( the single particle wave- function of 
h+ for the band e s n k in this periodic potential V(y) has 
the form 9^(x,y) = e i(k - x+k ^hi s nkv (y) with u s nky (y) = 
u n k (y + d) where k y is confined in the first Brillouin 
zone (BZ) k y £ [— Q/2, +Q/2], n is an integer, and the 
superscript s(= c, v) denotes the conduction (c) and va- 
lence (v) band. The wave function can be expressed 

as u n,kXv) = T,m C n,k v ( m ) eimQv i with the coefficients 
C n k determined by the central equation 



[h + (k x ,mQ- 



■k y )- 



]Ci k (m)+YViCt k (m-l) = Q 



i 



(3) 

with the condition C* fc (m) = S n ,mC>n,k ( n ) f° r ^ = 0- 
In this way, the electron band is folded into many mini- 
bands confined in the first BZ as illustrated schematically 
in Figs. 1(a) and (b), where the labels (n, s) represent 
different mini-bands. 

In the presence of V, the mini-band structure can be 
obtained numerically by solving the central equation. It 
is found that when the potential V increases the conduc- 
tion mini-bands shift down and the valence mini-bands 
shift up, eventually the two bands cross and form a in- 
verted band gap as shown in Fig. 1(c). Anticrossing be- 
tween the two bands occurs, which is induced by the in- 
teraction term A. A simple perturbative viewpoint helps 
us to apprehend why the periodic potential V tends to 
induce the band inversion. Take the band (0, c) for our 
illustration. For a weak V, the calculation up to the sec- 
ond order perturbation gives the energy correction for 
the band (0, c) 



<, k „\V(r)\K 



m,s,k' 



(4) 



c 0,fc 



here the summation excludes the band (0, c) itself, u s r 



and t s n k are the unperturbated eigen state and eigen en- 
ergy of the mini-band in the first BZ. Using the Fourier 

series of the potential, we have (w§. fc<i |V(r)| u s mk , 



V- m S 



U 0,k v \ U m,k l 



, °, k v 

V(r) couples the band (0, c) 

and other conduction bands (n, s — c) stronger than 
the valence bands since they come from same origi- 
nal bands and share the same spin indices. We have 



> 



O.fe' 



Thus we draw the 



conclusion that the perturbative correction is always neg- 
ative, Aep k < 0, whose effect is shifting the band (0, c) 
downward. Similarly, the valence band (0, v) shifts up- 
ward in the field. At a certain value of V, the two bands 
cross and open a gap again, i.e. a negative one. Simi- 
lar argument can be applied to other mini-bands. With 
increasing V, more and more mini-bands will cross near 
the zero energy point. Thus the band inversion occurs 
consecutively because of the shifting of the mini-bands 
from the conduction and valence bands. 

The formation of the inverted band gap is of topolog- 
ical and experimental interests. To gain a quantitative 
insight of this physical picture, we discretize the Hamil- 
tonian h + (k) in Eq. (2) on a square lattice. Without 
losing generality, we make the parameters material inde- 
pendent by setting A = —B = 1, D = and the lattice 
space a = 1. Eigen energies and eigen states are solved 
numerically for the periodic system. Following Thouless 
and his co-workers, [18| we come to calculate the Hall con- 
ductance of h+ in the band gap, which is equivalent to 
the Chern number of the filled valence bands. The Berry 
curvature for each band is defined as 1 19] 



dui 



dk x 



dui 



dk v 



Out 



dk v 



dui 



dk x 



(5) 

The first Chern number in the band gap is then com- 
puted by summing up Berry curvatures of the occupied 
bands and integrate over the first BZ C — ^ Yl n II ' 
^n,v(k)- The Hall conductance for h + is related to the 
first Chern number by a xy — ^-C.[3, [2(| In the ab- 
sence of the potential, the Chern number was found to 
be C = -[sgn(A) + sgn{B)}/2.^ The signs of the two 
model parameters in the term (A/2 — Bk 2 )a z of h + de- 
termine the value of C, or ±1, and further whether the 
system is topologically trivial or not. In the presence of 
V, it is found that the Chern number changes by 1 when 
the valence and conduction bands cross. By continually 
changing the potential period d and potential magnitude 
Vq, we obtain a phase diagram for the Chern number in 
Fig. 2 with the numbers in the boxes indicating the val- 
ues. We see that in this system the Chern number is a 
function of the potential period d and amplitude Vq. 

Examination of the evolution of the band structure and 
Berry curvature at the transition courses further confirms 
our claim. In Fig. 3 we show the first two transitions 
from C — to 1, and from C = 1 to 2. By increas- 
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FIG. 2: (Color online) The phase diagram of Chern number 
in (Vb, d) plane for the model with the Fermi energy Ef — 0, 
(a) A/2 = 0.01; (b) A/2 = -0.01;. The dashed line in (a) is 
obtained by solving Eq.(2) in the truncation approximation 
for m = and ±1 in C, s n fcy . Other parameters are A — —B — 
1, D = and the lattice space a — 1. 



ing Vb, the band gap closes and re-opens in Fig. 3(a) 
and 3(b). Correspondingly, the Berry curvature peak 
become sharp, and reverses its value dramatically, which 
accounts for the change in the Chern Number. At the cor- 
responding k points, the eigen wave functions also switch 
their band indexes quickly. By further increasing Vb, the 
peak of the Berry curvature splits into two sub-peaks, 
and the third peak with an opposite sign grows up while 
the Chern number remain to be C = 1. In the second 
transition, the third peaks reverses again just like the first 
transition, and the Chern number changes from C = 1 
to 2. Therefore accompanied by every "closing and re- 
opening" of the band gap the Chern number changes by 
one. 



According to the edge-bulk correspondence 21, 22], the 
first Chern number determines number of the edge states 
when the system has an open boundary. We take a strip 
geometry that is parallel to the x direction with open 
boundary condition along the y direction. We start from 
a topological trivial case of A > and B < 0, and scru- 
tinize its changing in the band structure when varying 
the periodic potential Vq. In the absence of the poten- 
tial V the system is insulating and topologically trivial 
indicated by a positive band gap. As Vq increases, the 
conduction bands shift down and valence bands shift up. 
Eventually the two bands cross, and a band gap re-opens. 
A pair of edge states appears connecting the valence and 
conduction bands as shown in Fig. 4b. Detailed anal- 
ysis shows that the wave functions of these states in- 
deed reside on the sample edges only. As Vq goes further 
up, more bands cross and more pairs of edge states are 
formed at distinct k points as shown in Figs. 4c and 4d. 
It is interesting to observe that the band gap at E = 
always retains a finite value at certain stages throughout 
the evolution. We can see the appearance of the edge 
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FIG. 3: (Color online) (a)-(d) Energy dispersion for the first 
valence and conduction band for various periodic potential 
strengths (a) V = 0.06; (b) 0.14; (c) 0.23; (d) 0.26. Color 
red indicates the dominant eigen state is spin-up, color blue 
indicates the dominant eigen state is spin-down. The other 
parameters are A/2 = 0.01, d = 30. (a) and (b) corresponds 
to the transition of C = —* 1. (c) and (d) corresponds to 
the transition of C = 1 — > 2. (e)-(h) Berry curvatures of 
the valence band around F point corresponding to the band 
structures (a)-(d), respectively. 



states corresponds to the change of the Chern number 
from C — to C — 3 by comparing Figs. 3 and 4. Thus 
modulation of edge state numbers is clearly reflected in 
the step wise changing of the Chern numbers. 

h+ and h- are the time reversal counterparts of each 
other. Each edge state m ^+ nas a couterpart 

|^_) = Q m h-i an d they form a pair of helical 
edge state. As a result, the large Chern number (C > 1) 
indicates multiple pairs of helical edge states in the sys- 
tem Ho in this periodic potential V. For a thin film fab- 
ricated on the substrate, the top-bottom symmetry will 
be broken due to the interface of the thin film and the 
substrate. This fact will remove the degeneracy of the 
spectra from h + and An off-diagonal term of struc- 
tural inversion asymmetry is added to the model Hq in 

(To 

<r 

up and down blocks (ctq is the 2x2 identity matrix). 



Eq. (l),[2j| AV = VsiA 



which connects the 
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FIG. 4: Band stucture evolution of a strip geometry in a 
periodic potential. Strip width 400, d = 20, A/2 = 0.03. (a) 
V = 0,C = 0; (b) Vo = 0.3, C = 1; (c) V = 0.6, C = 2; 
(d) Vo = 1.0, (7 = 3. Asymmetric dispersion for edge state 
originates from the symmetric breaking of the potential on 
the stripe with open boundary condition. 

Similar physics happens if the sample breaks the bulk in- 
version symmetry in quantum well. [24j ] In the presence of 
VgxAi numerical calculation still demonstrates existence 
of the edge states in the four-band model, which is char- 
acteristic of QSH phase. Thus this term does not destroy 
the QSH phase explicitly. 

Feasibility of experimental realization of this phe- 
nomenon depends on spin decoherence length in the sam- 
ple, and the fabrication of the periodic potential. The 
spin decoherence length was estimated to be 1 ~ 2/j,m in 
the HgTe/CdTe quantum well,0Jl]| and expected to be 
longer in a thin film of Bi2Se3.[23[ On the periodicity of 
the modulated electric field, the sub-lOOnm period of SSL 
fabrications was reported a decade ago, [25| and the 3nm 
periodicity for a graphene on Ru(lll) recently. [26] For a 
set of realistic parameters for HgTe/CdTe quantum well 
with a normal band gap A = 4meV,0] the magnitude 
of the potential with periodicity d — 50nm is calculated 
to Vq — 40meV for the transition of C = to 1 and 
Vo = HOmeV for the transition of C — 1 to 2. We specu- 
late that the modern techniques of superlattice make this 
electric field modulation possible. 

In summary, the band structure of a thin film or quan- 
tum well are folded into the mini-bands in the reduced 
BZ by a periodic potential, and can be modulated such 
that the conduction bands shift downward and the va- 
lence bands shift upward. Each process of the band gap 
closing and re-opening will change the Chern number by 
one. As a result the number of the helical edge states 
will increase or decrease by one. This demonstrates the 
possibility of the electric field modulation of topological 
orders in the thin film semiconductors, which is analogous 
to the integer quantum Hall effect in a strong magnetic 



field. Direction of electrons of the helical edge states is of 
interests in quantum information and quantum process- 
ing. Controllable number of the helical edge states will 
pave an alternative route for application of edge state 
physics in the future. 
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